VIBRATII 
NELINIARE 


VIBRATIILE LIBERE ALE UNEI BARE TIMOSHENKO 


IN CONSOLA 


1. INTRODUCERE 


Ecuatia vibratiilor libere transversale unidimensionale este bazata pe ecuatia clasica 
Bernoulli-Euler care este inadecvata pentru modurile inalte de vibratie. Problema vibratiilor 
transversale ale unei bare Timoshenko se afla in atentia multor cercetatori care considera 


diferite metode de rezolvare. Bruch si Mitchell [4] considera o bara Timoshenko cu o masa 
aditionala la un capat. Abbas si Irretier [1] folosesc metode de tip Runge-Kutta-Fehlberg 
pentru determinarea pulsatiilor precum si metode experimentale. Gutierrez s.a. [5] folosesc 
metode de tip Rayleigh-Ritz optimizate. Benerjee si Williams [3] se ocupa de vibratiile 
transversale in spatiu ale unei bare Timoshenko cu ajutorul matricei de rigiditate. Abramovich 

2] se ocupa de vibratiile unei bare compozite folosind metode exacte de integrare. Yuan si 
Dickinson [8] folosesc arcuri artificiale in studiul vibratiilor unui sistem de bare. 

In cele ce urmeaza se obtin si se rezolva doua ecuatii diferentiale in doua variabile. Se 
folosesc conditii la limita pentru bara in consola, se obtine ecuatia frecventelor si cea a 
modurilor normale. 

Se considera un exemplu numeric din care se deduce ca inertia de rotatie si deformatia 


de rasucire nu se pot neglija intotdeauna. 
2. ECUATII DE MISCARE, CONDITII LA LIMITA, SOLUTII 


Folosind principiul lui Hamilton: 


ГЕ (Е, – u)de = 0 (1) 


“a 
unde ғ; si £, sunt doua valori fixate ale timpului, E. energia cinetica a barei Timoshenko, U 


energia sa potentiala, putem obtine ecuatiile de miscare. 


Energia cinetica este data de relatia: 
_ 1н /(Фү 1 ri әү 
Е, ЕСЕ РА x) dx +>, et (=) dx (2) 
iar energia potentiala de formula: 


u - if ot (2) de-i fi kac(Z — y) ах G) 


unde y(x,t) reprezinta deplasarea transversala, u(x,t) unghiul de incovoiere, A aria sectiunii 


transversale, p densitatea de lungime, I momentul de inertie geometric al ariei sectiunii, E 


modulul de elasticitate, G modulul de rigiditate, k un factor numeric de forma al sectiunii 


transversale, x€ [0, L] este abscisa unui punct curent de pe bara de lungime L. 


Din principiul lui Hamilton (1) se obtin ecuatiile diferentiale: 


ay 
as 


au ay = 
E73 + kac (2: —ар)-р17-# = 0 


с [2х му дэ. 
ке(22-2)-ғА -о 


Folosind notatiile а/а:=*, 8/8х-” si 


ecuatiile (4) si (5) devin respectiv: 
a^i" 4 = (y' 2- 4) = a! bu? -0 


y'-y'-ay- 0 


(4) 


(5) 


(6) 


(7) 
(8) 


Eliminand + si y intre ecuatiile (8) si (9) se obtin doua ecuatii diferentiale in y si 30: 


ӘУ (a? +b?) ЗУ e un zx = 
(a? + b?) + apt 22 =0 


ax àx?art r art 
д*%ф 2 2, 8% BG, 11289 _ 
ae Eee эр tab аш? 


Conditiile la limita impuse pentru bara in consola sunt: 
у(0,0- (0,2) = 0, w(Ltü)-y'(Lt), w'(Lt)=0 
Aplicand metoda separarii variabilelor ecuatiilor (9) si (10): 
у(х,г) = Х.(х)Т7(ғ%), Ч(х,т)-Х.(х)т(ғ) 
Se obtine: 


Klx) = €C,cosqx + C, sin qx + C4coshsx + C, sinh sx 
‚ z x іл, 
Х.(х)- AE (C, cos qx — C, sin qx) + TE (Cs cosh sx + C, sinh sx) 


unde 
4 TP 21 IE өрен 4 
2g^—p HG + b*)* + 42] + p*(a* +b“) 


= ть = = 1/2 = т 
2s? =p [>°(а +?) + 45| —p*(a* + b°) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


(15) 


(16) 


T(t) + p*T(t) = 0 (18) 
si conditia b >> арт pentru cazurile practice. 
3. ECUATIA PULSATIILOR PROPRII SI MODURI NORMALE 


Conditiile la limita (11) in X, si Х- sunt: 
Х,(0) -Х.(0) 20; ХІ(І)-Х.,(1); x;(z)= 0 (19) 


Ecuatia pulsatiilor proprii devine: 


c 458, =) cosqL cosh sL ++ 


хээг Бич гн - sin gLsinhsL — 2 = 0 (20) 


Radacinile FP., i= 1,2,3,... ale ecuatiei (20) dau valorile proprii ale problemei 


ы» 
considerate. Modurile normale corespunzatoare X;, si X; sunt de forma: 


= | 7 zi à ыг apja 
Х.-К,; (соз q;x t Ө,сіпд;х ор coss.x + =L E = B,sinhs,x ) (21) 
x, m K sing x — B,sinq;x —— ті sinh s, x+ Өтлтл 18, coshs,x | (22) 
ai е : а2р?- 
ипде: 
в = cosh s;L—cog q;L (23) 
i кіп Z; тыч нр 
qita p-ar) 
Cu X4; si Х-; cunoscute, deplasarea transversala si unghiul de incovoiere pot fi scrise 
sub forma: | 
¥(x, t) = 2; а,Х, е Рр) (24) 
х, t) = Ej b, Хе етер (25) 


unde constantele а, bj, e; si В, se determina din conditiile initiale. 


4. EXEMPLE NUMERICE 


In multe lucrari se afirma ca neglijand inertia de rotatie si deformatia de rasucire, 
pulsatiile astfel obtinute nu difera prea mult de cele obtinute prin luarea lor in considerare. 
Dar pentru unele bare (de profil ,I’, ,T’, ,U',...), aceste frecvente difera considerabil. Pentru a 


justifica acest lucru, presupunem о bara de otel de profil „Г. Presupunem 


Е=2,1.107[М/ст°]; G-3E/8; р-78-10 ?[Kg/cm?] А-30 [cem]; 


k=1/4,4. Se obtin а? = 4,27 - 107" [(s/cm)?]; b? = 3,64. 10 1 [(s/cm?)?]. 


I-625[cm?]; 


Primele patru pulsatii pentru diferite valori ale parametrilor L si r^ corespunzatoare 


barei Timoshenko (T) si Bernoulli-Euler (B-E) sunt date in tabelele 1-8: 


jP 100 200 300 
r 
125/4 329,616 83,41 39,17 
125/6 289,86 68,48 29,54 
125/8 234,24 58,56 25,58 


Tabelul 1- Pulsatia p, pentru bara B-E 


100 200 300 

2 
125/4 295,17 79,58 35,81 
125/6 248,86 65,46 29,34 
125/8 220,29 56,90 25,45 


Tabelul 2- Pulsatia p, pentru bara T 


L 100 200 300 
r` 
125/4 2040,59 510,14 226,73 
125/6 1666,13 416,53 185,12 
125/8 1442,91 360,72 160,32 


Tabelul 3- Pulsatia p. pentru Бага B-E 


100 200 300 
p 
125/4 1352,84 443,66 212,04 
125/6 1223,36 377,90 176,86 
125/8 1125,81 334,75 154,87 


Tabelul 4- Pulsatia p. pentru bara T 


L 100 200 300 
p 
125/4 5713,71 1428,42 634,85 
125/6 4665,22 1166,30 518,35 
125/8 4040,20 1010,05 448,91 


Tabelul 5- Pulsatia р; pentru bara B-E 


100 200 300 
p 
125/4 2933,17 1081,17 549,06 
125/6 2735,75 951,66 468,21 
125/8 2580,00 860,31 415,08 


Tabelul 6- Pulsatia p, pentru bara T 


L 100 200 300 
г? 
125/4 11196,60 | 2799,14 | 1244,06 
125/6 9141,98 2285,49 | 1015,77 
125/8 7917,18 1979,29 | 879,68 


Tabelul 7- Pulsatia p, pentru bara B-E 


L 100 200 300 
кі 
125/4 4529,18 1825,70 978,83 
125/6 4330,40 1652,39 854,36 
125/8 4158,75 | 1522,38 768,15 


Tabelul 8- Pulsatia p, pentru bara Т 


Se observa ca efectele inertiei de rotatie si deformatia de rasucire nu pot fi neglijate 
intotdeauna. 
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VIBRATIILE LIBERE ALE UNEI BARE CU UN NUMAR ARBITRAR DE MASE 
CONCENTRATE 


1. ECUATII DE MISCARE SI SOLUTII 


Utilizand distributia Dirac, legile dinamicii, ale fortelor elementare si ecuatiile 
generale de elasticitate, vom obţine următoarele ecuaţii ale vibratiilor libere de incovoiere a 


grinzii finite, uniforme. 


= — [oA - Ez C m,O(x— х0123-0 (1) 
ам 

YT- pe 55-0 (2) 

М = —EI E (3) 

T= kaG(2 -а) (4) 


in care y(x,t) reprezinta abaterea transversala, «(х,Ё) unghiul de incovoiere, T(x,t) forta 


transversala de forfecare, M(x,t) momentul de incovoiere, E modulul de elasticitate, С 
modulul de rigiditate, Т momentul de inertie geometric al ariei sectiunii, A aria sectiunii 


transversale, p densitatea de lungime, k un factor numeric de forma al sectiunii transversale, 
m. (i = I,m) o masa care actioneaza in punctul de abscisa x,¢ (0,1), | lungimea grinzii. 


Eliminand parametrii M si T intre ecuatiile (1), (2), (3) si (4) obtinem: 


EIES + кас (Ë - a) - pte — 0 (9) 


Prin derivare si eliminand parametrul & din ecuatiile (5) si (6), se obtin urmatoarele 


ecuatii diferentiale in y si а: 


Et[oA* E... түйіс sal 25 8 п, m,ó(x _ ч) 2 


EI {р + КАФ 
- ER 1m, 6" (x — x, ДЕ -- шім» 25 mjb(x-— x, DE 


- =* D. NE n Plat". 


кА dz аз 


(7) 


ЕгА+ Т. түё(х ah ete 


= {рї + = -+ [pA + Үл т,б х,)] 2° 


х®ё@г® L= 


F at a 
+ 7 [pA + Dy m  х- ха = 


= – Ўн m(x-— x)22 
(8) 


unde 2 (х- x,) si "(x — x,) sunt prima, respectiv a doua derivata a distributiei Dirac. 
Metoda separarii variabilelor poate fi aplicat in ecuatiile (5), (6), (7) si (8). Astfel, 


separam variabilele in ipoteza ca: 


y(x,£) = X, (x)e'** (9) 
a(x,t) = X;(x)e'*, i= 4—1 (10) 
Folosind notatiile 

а. 825 ж.т Q2. BAM 2298 

a ах ' ч к” b gt E (1) 
relatiile (5), (6), (7) si (8) devin: 

а? Ху + (a2c? — b)X, + b? X? = 0 (12) 
X; + aX, += yn =, m,0(x — x,)X, — X; = 0 (13) 


27 + (a? + c2)X + (a?c? — b2)X, 


а2с2-в° 


+. mire +2 zy fa, ме, [d(x — x,)X,]" = (14) 


| . шоо 25 z z_ = ел = I 2 
XD" + (а? + саулы + (a?c? = БЭХ, + Е = Li 09 т,б(х х,)Х + 


a n T B n Г 
E mx —x,)X; = "OLI (x — x;)X, (15) 


Aplicand o transformare Laplace ecuatiilor (14) si (15) se obtine o ecuatie elgebrica de 
forma: 


[gf + (a?  e?)q? + а?е? — 57]X,(q) — q*X, (0) — q*Xi (0) — ах! (0) — x;" (0) - 


—(a? + c?) [ax, (0) + x; (0)] + ——— 254 Ун зайн 
ізі 


жау me Ua X0 = 0 (16) 


[q*+ (а? + c2)q? + а?е? — b2]X;(q)— ala? + а? + с?)Х,(0) – (а?+а? + c?)x;(0) — 


—-qX; (0) — xz (0) j s іі тп; e 9i OG (x; )+ 
утте о) = em see [aX (к) — ЕЗІ (17) 


unde Х, (9) este imaginea lui X, (x) prin transformarea Laplace (k=1, 2). 


Cu ajutorul expresiilor (6) si (17), transformarea inversa ne da: 


x,(z) = TIE G cosh ex + a? cos [jix)X,(0) + (E sinh ах + Tm ex) X,(0) + 
+(coshax — cos Bx)X4 (0) = Č sinh ax — 2 sin х)Х( (0) + 
y X (x) Ë diae = sin B(x - š 2) нбс) 
= 
_ «уу, m,X, (x, = 2 sinha(x — x,) +Ë — sin B(x — х.)Н(х x) | (18) 
X,(x) = = G coshax + a? cos B£x)X,(0) + e sinh ах + Tsin 2 Х:(0)-- 


1 
+(cosh ax — соз f x)X7 (0) + sinh ах- z sin gx)X; (0) + 


[ас 5 х. a(x — x,) - sin B(x = xD) шін 


іші 


Е иб m;X7 (x) [sinh a(x — x,)— ;sin В(х- «ӘН (ж —x,) + 


+ "ys. m,X4(x,)(cosh a(x — x,) — cos B(x — x;)) — 
= ух. m Xi (x) C sinh a(x —x,) -узш B(x - ЛДЕ — xj) (19) 


unde 
242--а — с? +4) (а? —c?)? + 462, 282 =a* Rc +4 (а? - с) + 482 (20) 
si Н(х- x,) este distributia Heaviside. 
Dar, doar jumatate din constantele X, (0), Х,(0), Х,(0),Х, (0), k = 1,2 sunt 
independente in relatiile (18) si (19). 


(a* 9 5*)x; (m) -a* x; (а) 


х,(0) = ‚ X;(0) = a^X,(0) + X7 (0) 


Bac (21) 
x: (0) = a°X1(0) +f" (0) + X7 (6) = (b? — a?c?)X, (6) — c2x7(0) 
2. CONDITII LA LIMITA 
Conditiile la limita necesare si suficiente sunt urmatoarele: 
X,(0)=X,(0)=0, (0) = X,(D, Х(1)-0 (22) 
Pentru acest tip de fascicul, din conditiile (22), obtinem: 
X (0) = У, Cy X,(x,) + У. C; Xy (x;) + 
+ Lica Са) + Dh Cy (ХХ: (xi) (23) 
X3 (0) = 2-4 CX, x) + Ty Ce Xi (ху) + 
ERG CaL) + EI Cg: (х,)Ху (zi) (24) 
unde : 
m;a? | = ñ 
С = prid coshal + a? cos Bl — ac? sinhal — Bc? sin fl +c2) - 
E E! 1 
-casha(i — x,) + (В? — с2) cos f (1-х,)- b° (—sinh al — 8 sin !)- 
- z sinh a(l— х.) + B sin B (L— x;)) (25) 
p № aC ынай рев 
^ = —— a2(— sinhal — = : — x, —x,)- 
2 pAAl `a B ; 


ü [sine -х)- әзер u= х0| 


- (В? coshal + a? cos Bl — ac? sinhal - Bc? sin BI) (26) 


E m,(b? — a?c?) 


pAA 
- (В? coshal + a? cos BL — ac? sinhal — Bc? sin BI) — 


Cs (ГЕ sinha(1— x,)— jim ға- х) j 


-a2€ sinh al — аш [cosha(1 — x,) — cos (1— х1) i: 


2 
т;а 


С. = 
ч pAA 


ec sinh al — езіп Bi)lcash a(! — x;) — cos Ë (1 — x,)— 


si ЕС a(l—x)- sin p(t = х). 


- (В? coshal + а? cos Bl — ac? sinh el — Bc? sin BI) (28) 


m; 
рАА 


С; = ((a* + b*)(coshal — cos #1)[(a? + с?) cosha(l — x,) + 


+(В? — с?) cosB(l — x,)] — ^ [(8? — a?) coshal + (а? + a?) cosh ВП. 
-a sinh a(l — x;) + B sin B (1 — x;)]) (29) 


NEL. 8 2 2.. 23 1 24,1 ВП. 
Сы; ақ” [Cg ) cosh ai + (a^ + c^) cos ВИ 
4. 
(созһа(1- х,)- соз (I — x,)] — — —(coshal — cos B D) - 


- [sinh a(t — x) — sin В(1- xj]! (30) 


m;(b? — a*c?) fat + bš 1. 1. 
бүс — 201 k - соғ 1) [sha Qi —x;) — газа": = | = 


—[(8* — a?) coshal + (a? + а?) cosh ВП. 
[coshe(1 — х) — cos (1-1) (31) 


аёт; 
Са = 
рАА 


{[(8* — a?) coshal + (a? + а?) cos 81]. 


аё z 
[cosh а(1- x,)— соз B (1 — x;)] — 


= (coshal — cos Ё Г) - 


- [sinn a(l— x;)— геп В (1 = » 


A- — cosh al sinh al — B° (В? — a*)(coshal)* + 
ES TA cash al sin BI + (ас — а? — 2b?) coshal cos fil + 
желе t sinh al cos Bl — a (а? - а? Icos gi)? + 


B*e (a +a" \-а*—һ° 


2 “кіп Bleos Bl 


Relatiile (21) relativ la (22) sunt reduse la: 


nr a* t bš t тг 
Хү'(0)= ---т-Х(0); Х0) = x7 (0); 


Xi(o)--5x(0 (0) = 2! (0) 


Solutiile (18) si (19), relativ la (23), (24) si (34) devin: 
X (xz) = Em, IU, X x) + Dy U, X (x;) + Dy 5,Х.(х,) + 
+ Ei ИХ: (xj) 


Ха (х) = Eia V, X, (x) + TE, Vy Xi (zi) + TE, Vy X; (x) + 
+ Dia Vi X; (x;) 


_ 1 p? a? + ы | a a? + b^ 1 
Ua ара la а sin ах + BP pa sin Bx Cy t 


+ (cash ar — cos f Ca + [b° sinh cado T с 
(coshax — cos x)C,; zA Gs a(x — x;) yt (x —x;}) 


-a (e sinh a(x — х;)- шан sin В(х — xg —x)] 
1 B a 
ы PFa sin ax + % sin х | C;, + (coshax — cos B x)C,; — 
= Ра sinhax — ¿sin х), | 


(32) 


(33) 


(34) 


(35) 


(36) 


(37) 


(38) 


1 2 2 
Шз; = "EY (£ sin ax + sing) C3; + (coshax — cos B x)C,; — 


- £I sinhax — ;sin Bx)C;,| (39) 
7 - (5 LT Je tos Bx)C 
қ sin ax sinfx | C, — cos B x)C,. — 
^i EE ny B 4; + (coshax — cos B x)C, 
= неме sinh ax — 3 sin Boc] (40) 
2. 2 21 „2 2 
у; = EY (En ax + E 2 sin Bx)C;; — = (созһах — созВх)С,: + 
+ = [cosha(x — х.) — cos Ё(х—х,)]Н(х— 2) (41) 
PER. 2 28 „2 
Va = = - gi (Ea ax + = 2 = imi Bx)G;; — (cosh ax — cos Bx)C;; — 
- — * sinh a(x —x,)— asin B(x — х) Н(х-х 2] (42) 
z a 21.2 2 
Va = = 2 gi (4 = E_sinh ax + = T sin Bx)C,; — "(cosh ax — cos #х}Су + 
D “|: sinh a(x — x,) — a sin B(x — x; Ә|н(х — x, 2] (43) 
2 _ 2 ЭГ Y 
Vy = = — (En ах + E - ы sin х)С,,- 2 (созһах - cas f x)C,, — 
-Tn [E sinh a(x - а) asin B(x —x)] nr — x;)} (44) 


Derivand in ecuatiile (34) si (36) obtinem: 


X (х) = Жы, Ur X, ұ)) + ET, Vat + Хе. Utal + 
+ iza Uşa (x) (45) 


X; (x) = X, Vit (x) + Уш (a) + Eta Vata) + 


7 


+ Dik Vu Xa (к) (46) 


Š 


Pentru x = x, (j = 1,n) in relatiile (35), (36), (45) si (46) si pentru solutiile netriviale 
Х,(х,). Х.б), Xz(x;) si XZ (x), (i = Tn) determinantul coeficientilor trebuie sa se 


anuleze, pentru care se obtine: 


(Ui: (х;) - бн) (U2; (x 9) (U3; (х;)) (0. (5) 
| 048) (0 08)-80 Wa) 05) 1., 
Qux) EU CRED — &) GAER) 
(АО) | 0100 (Z) 0 (5) 8), 


cu í = 1,m si j = 1, п. 
Fiind stiute X, (x). X. (z), putem determina abaterea totala, unghiul de incovoiere si 


momentele. 


3. CONCLUZII 
In lucrarea de fata sunt evidentiate urmatoarele: 
(a) solutiile sunt obtinute din doua ecuatii diferentiale pentru o incovoiere totala, 
respectiv pentru unghiul de incovoiere, 
(b) constantele din aceste solutii sunt legate de una din cele două ecuaţii originale din 
care cele două ecuaţii diferenţiale, de mai sus, sunt derivate, 
(с) condiţiile la limită prescrise sunt omogene. 


Astfel, se obtine frecvenţa şi ecuaţiile normale pentru modul fasciculului de consola. 


VIBRATII TRANSVERSALE ALE UNEI GRINZII IN CONSOLA CU MASA ГА САРАТ 
SUPUSA EXCITATIEI ARMONICE LA BAZA 


1.Introducere 


Problema determinarii influentei unei mase de varf asupra comportamentului de vibratii transversale a unei 
grinzi ancorate la celalalt capat a fost studiata de mai multi autori. Masa de varf,care este utilizată în multe 
desene sau modele structurale, ar creste flexibilitatea sistemului si poate modifica comportamentul vibratiei 
sistemului. 

Hoppmann a studiat o vibratie speciala fortata a unei grinzi care transporta o masa la mijloc, fara a investiga 
problema valorii initiale generale.Metoda pe care Timoshenko a folosit-o pentru a rezolva un caz special al unei 
probleme a vibratiei longitudinale a unei tije cu o masa atasata la capatul ei, nu este usor de aplicat intr-un caz 
mai general al vibratiilor transversale ale grinzilor. Chen a investigat problema dintr-un punc de vedere 
matematic. El a folosit metoda separarii variabilelor si a determinat functiile proprii ale ecuatiei diferentiale 
ordinare rezultata, dar nu a luat in considerare nici nelinearitatile sistemului nici excitatiile de la baza.Grinzile 
subtiri cu o masa de varf si excitatii de baza sunt reprezentate de un model matematic simplu, ca si majoritatea 
structurilor antenelor, fiind investigate de To, care a evaluat frecventele naturale si modul lor corespunzator 
pentru modelul linear.Laura a studiat frecventele naturale ale unei grinzi prismatice fara cleme, avand o masa 
concentrata la capatul sau. 

Este cunoscut faptcul ca multe cercetari au investigat problema grinda-masa, cautand o solutie armonica pentru 
ea, dar cele mai multe au utilizat ecuatii liniare diferentiale cu prezumtia ca exista un raspuns armonic. 

Este evident ca modelele liniare pentru aceasta problema ar avea un raspuns periodic, fie o forta armonica fie o 
excitare armonica a bazei.Acesta nu este cazul pentru modelele non-liniare.Prin urmare, trebuie dovedit faptul ca 
ecuatia diferentiala non-liniara a sistemului are un raspuns periodic, iar conditiile periodicitatii trebuie sa fie 
dezvoltate de asemenea. 

In cele ce urmeaza, un sistem interesant de modelare a unei grinzi in consola simpla supusa unei miscari 
armonice la sprijinul sau este de analizat. Ecuatia non-lineara de miscare a fost derivata intr-un mod general, iar o 
abordare matematica a problemei este prezentata.Unul din cei mai importanti pasi in teoria oscilatiilor non-liniare 
este de a stabili un criteriu suficient pentru a garanta existenta unor solutii periodice.Cele patru metode cele mai 
frecvent utilizate pentru a dovedi existenta unei solutii periodice sunt dupa cum urmeaza: 

а) metoda continuitatii analitice 

b) procesul de echivalare al coeficientului Fourier 
C) aplicarea de teoreme de punct fix 

d) metoda seriilor de puteri 

Pentru rezolvarea ecuatiei diferentiale de miscare am folosit metoda homotopica, iar rezultatul este comparat cu 
solutia numeric rezolvata. 


2.Ecuatia de miscare 


Figura 1 reprezinta o bara subtire, cu lungimea /, si masa de capat m, studiata in cadrul de referinta cu 
coordonatele x-y. Masa grinzii se presupune a fi mica si infima comparata cu masa concentrata si, prin urmare, 
grinda este considerata ca un element de incovoiere.Suportul grinzii este supus unei miscari armonice v, in timp 
ce grinda este limitata sa oscileze numai in planul x-y. 

Pentru derivarea ecuatiei miscarii pentru sistem, dezvoltam expresia Lagrangiana.Energia cinetica, Ec a 
sistemului este datorata masei de varf.Prin urmare, energia cinetica a sistemului este 


Е = my2 +(Х + „?| 
€ 12 

Unde Y este deplasarea laterala a lui m, X fiind deplasarea axiala a lui m, si v este deplasarea prin excitare fata 
de sprijin. Nu exista nici o energie cinetica din cauza inertiei de rotatie pentru punctul de masa. 

Energia potentiala totala, Ep consta in energia de deformare cauzata de indoirea grinzii, si energia potentiala 
gravitationala a masei de varf. (2) 


Putem pune in relatie deplasarea laterala y cu centrul de masa, Y=y(I,t), folosind urmatoarea relatie bazata ре 
primul mod al formei liniare o(x) : 


Unde 
Ф (x) = 1 - cos LE = D| (5) 


Inlocuirea ecuatiei (4) in ecuatiile (2) si (3) rezulta: 


| I y ! (6) 
Е = o | а? т dx - mgX BEN tA - тех 
p 3 |, 2 
5 i (7) 
i 2 
1 
X = — У 2 4. ) dY = ВҮ 2 
2 | d J 2 
0 
Unde 
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Dupa се integram ecuactia (7), obtinem: 
X-BYY (8) 
Si substituind in ecuatia (1), L din Lagrangian va arata astfel: 


(в? 22 ) 1 2 1 2 
L=E -57 Boy Y? + ү? + v” = 2BYYv - ВУ А + >-mgBY (9) 


Aplicand ecuatia Іші Lagrange se obtine urmatoarea ecuatie diferentiala de miscare pentru sistem: 


EI 
( B үд) х В а EZ - gB - Bv(t) ҮҮ -0 (10) 
m 


) 


Aceasta este o ecuatie ordinara parametrica non-liniara, care, prin neglijarea termenilor non-liniari, se reduce la o 
ecuatie diferentiala de tip Mathieu.Cu toate acestea, daca luam in considerare doar excitatia armonica simpla a 
bazei de forma: 


v = Qcosat (11) 


Apoi urmatorii parametrii pentru ecuatia (10) sunt generati: 


T Ё 1 п 7 1 
= Б = Dž] 5 B- Б = о E v = -Qo?.cos ot (12) 


Introducand unrmatoarele simboluri si variabile fara dimenisuni 


" ВҮ !- du 
u=BY ; u-BY/lo; "7 EN саа dq 21-6 (13) 
© 

putem reduce еспайа (10) la forma: 

и" + Cu + wu" + aha: + Du:cos = 0 (14) 
unde: 
E ERN ESSA ЗЭЭ 

р-ОВ , =з 032) ; u(0)=B ;u(0)-0 (15) 

m 


[0] 


3.Solutia periodica folosind metoda homotopica 


Noi construim o familie uniparametrica de ecuatii 


ТТ ТОТЫ ТТІ 


Cu conditiile initiale 

dU(n. p) (17) 

U(n,p) = В = =0 pentru п -0 

Unde pe[0,1] este ип parametru , U(n,p) “este acuma o functia de n si p pentru ca О(п,р)=ио(п). 
La p=1 ecuatiile (16) si (17) sunt asemanatoare cu ecuatiile (14) si (15), pentru ca U(n,1)=u(n). Ca parametru р 
variaza intre Osi 1, U(n,p) si A(p) au formele: 


U(n.p) = u (n) + pu | (п) + рга, (п) +... (18) 
1 2) 2 
А(р) = С^+ (4 + (4 Ги” (19) 


Unde functia uk(n) si constanta 19 sunt functii nedeterminate 
Fiind p=1 si notand О(п,1)=и(п) si A(1)=A se obtine: 


u(n) =u (n) +u (1) + u, (m) +... (20) 
= С^ + pur" эс (21) 
Inlocuind ecuatiile (18) si (19) in ecuatia (16), rezulta: 
u (m) + Cu (m) = 0 (22) 
u(r) + Cu (m) = fug) (o) — ge (u^) = Рцту-совп = (a) Put (23) 


uy (a) + Cur) = Aug) чү (т) —2u (nu (mu (m) — 2u ТТ rit) — (Uu 01) — Du (n)eosn - (Т ШЕТПЕ ТЕТІН 


(24) 
Presupunand ca aproximarea initaiala a ecuatiei (22) si tinand cont de ecuatia (17) obtinem: 


u' (m) = B-cos(Cn) (25) 
Substituind ecuatia (25) in ecuatia (23) avem: 
Д | Д 
шб) + Cup = БХ _p2)_ ТЯ22 4 Ley + p2).cos(3on) - + BD(cos(C+ Dn + cos(C - 1) n) 


(26) 
Pentru a evita prezenta unui termen secular trebuie : 


1 1 2 2 
(4) Fs "us (3 -p ) (27) 
Obtinem urmatoarea solutia generala a ecuatiei (26): 


no m m 21225 1 anf cos(C+ In со8(С-11)) (28) 
41(0) = А:соѕ (Ст) + B:sin(Cn) 53 81 +В ) cos (3-Cm) + 2 zi TES 2621 


Unde constantele A si B vor fi determinate din conditiile u1(0)=0, u'1(0)=0 si obtinem: 


a = (+ g?) +e ; B-0. (29) 


4C -1 
Substituind (29) in (28) obtinem: 


2 2(2C + 1) 2(2C - 1) 


u (t) = 1p * 82) cos (Cn) : —cos (3:Cn) + BD 
32 4C -1 


cos (Cn) x cos (C + 1)n _ cos (C - 1) | 


Drept urmare avem solutia aproximta de ordinul intai a ecuatiei (14): 

C С+1 С-1 
и (t) = B:cos(Cn) + В [| + p2).cos(Cn) = cos(3-Cn) E pD cos(Cn) + cos(C + Dn SE, cos(C - Dn 
1 32 ас? =f 2(2C + 1) 2(2C - 1) 


(31) 


Rezultatul este comparat cu solutia numerica rezolvata in [7] folosind urmatoarele caracteristici: C2 = 1,5; D = 
0,05 si D = 2, conditiile initiale: U (0) = 0,5, и '(0 ) = 0. Figura 2 ilustreazá istoria de timp a sistemului de excitatie 
vertical. 


u(t) 
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p 
1 
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n 
Fig.2 Evolutia in timp pentru deplasarea laterala a masei de varf si excitatie verticale 


Pentru a demonstra validitatea conditiilor periodice adaptate, ecuatiile diferentiale de excitatie verticala au fost 
rezolvate numeric. Figura 3 arata ca sistemul excitat vertical are un raspuns periodic, care indeplineste conditia 
de stabilitate [7]: 2 D 

тахи) <1--- 


c 


Totusi, pentru C2 = 1,5 si D = 2, in afara regiunii de stabilitate ar ауеа un comportament instabil. Pentru valori in 
afara intervalului de stabilitate C2 - 0,5 si respectiv D - 2, performanta este foarte instabila. 


a 
Fig.3 Caracteristicile planului fazelor pentru excitatii verticale. 


4.Concluzii 


Vibratia transversala a unei grinzi cu consola care prezinta excitatii armonice la baza este caracterizata prin 
ecuatia diferentiala ordinara non-liniara parametrica. Comportamentul grinzii este fortat de amplitudinea dinamica 


si frecventa excitarii la baza. Limita superioara, care garanteaza existenta unui raspuns periodic, poate fi gasita la 
expresia de inegalitate (32) pentru o valoare data a perioadei de excitare. 


VIBRATII LIBERE АГЕ UNU SISTEM DE BARE CURBE 
IN TREPTE 


1.Introducere 


Metoda arcurilor artificiale este propusa in cadrul metodei Rayleigh-Ritz pentru 
studiul pulsatiilor proprii. Alegerea functiilor de coordonare se poate simplifica 
considerand problema data ca un caz "limita" a unei probleme la limita: rigiditatea 
arcurilor artificiale se allege suficient de mare ca sa aproximeze sustinerea rigida. Astfel 
de probleme au fost studiate de Kato, Mizusawa, Filipich si altii, Warburton si Edney , 
Gorman, Yuan si Dikinson . Este important sa se puna conditia de continuitate pentru 
deflectie si unghil de inclinare la jonctiunea dintre doua bare vecine. Functiile admisibile 
se aleg astfel ca sa satisfaca conditia la jonctiune si la limita. Rigiditatea arcurilor 
artificiale se alege mare, astfel incat se aproximeaza rigiditatea sistemului de bare. Se dau 
exemple numerice care se compara cu solutia exacta . 


2.Descrierea metodei 


Sistemul de bare este format astefel incat prima bara este incastrata la un capat si 
conectata rigid cu a doua bara. Raza de curbura a barei I se noteaza R, si subaintainde 
unghiul pi . Tangentele in punctele de jonctiune intre elementele adiacente sunt aceleasi. 
Se vor considera numai vibratiile in planul sistemului de bare. Deplasarea tangentiala a 
barei curbe i se defineste prin variabila v, iar fiecare bara se considera omogena. 


Functiile admisibile care satisfac aceste conditii trebuie alese pentru fiecare bara 
componenta. 


v, (a) =>) ав” i=1..p а, =a, =0 (1) 


Energiile potentiala si cinetica maxime pentru fiecare bara sunt summate si obtinem 
respective 


14 В, [v'(a) v (o), | 
ээ үеээ a 00: 
V ax po XE R° (o) a (2) 


Žo "Sala a ІМ (a) ЛЕЦПС У] Ча (3) 


unde p, este densitatea barei I, E, este modulul lui Young de elasticitate, A, este aria 
sectiunii transversale , I, este momentul de inertie de arie al barei fata de axa neutral, © 
este pulsatia naturala a vibratiei sistemului considerat. 

Vom introduce trei tipuri de arcuri artificiale: unul de translatie normala cu rigiditate 
Ку, unul de translatie tangentiala cu rigiditate К, si ultimul arc de rotatie cu rigiditate К». 
Energia de deformatie din arcuri este : 


њ 25 [v (B)- Vial 0] +5 ЭЭ [w(B)-v.( 0] +1 S DORIA учым ш 


si energia totala potentiala este V,,, = Vu +Vp- 


Principiul de minimitate Rayleigh-Ritz este folosit pentru rezolvarea problemei de 
valori proprii. In cazul nostrum putem scrie : 


i=1.p  j=1.m (5) 


Substituind V,,, si T a in relatia (5), se obtine: 
0 Р ! 3 ñ 3 " 
Ga. > Ch, B; | UR (о)[и a.) + v, (a a) | J eR ЕЛ p.) eat ЛЕ 
ij К=1 0 


Sh | v, ( (В) = wa | (0)| ЭЭГ [ (с.)(у, (В, ) +V, (В, )) = (6) 


Вы (0)(у. (0) еу, (0)) | суа bh (Ое да = 0 


k=l Pk 0 U, (а) 


unde am folosit notatiile: 


, EI ‚ РА P L. 
= =; L,=|R(a)da ; U,(a)=—— ; 
i EL с pA | ( ) ( B.R (а) 
k = Кы, k = к 9 k. = kL, О? pA i 2 
М EL? “ТЕ” * EL 7 ЕІ 


In cele ce urmeaza vom considera R, = constant si convenim : 


kyo = Kro = Kao = К -k 


Tp 


-К,-0, К-і.т (7) 


Мр 


In acest са?, ecuatiile (6) sunt de forma : 


У | (4, т O Bá ) qi, + са | + Daia + О»а > E Рз таз = 0 (8) 


р 
К- 


Ын 


i=l.n; j-1.m 


Am notat : 


д, =ош о-ва (30-20-3403; 


(1-1/6-1),, 


ЕЗ” E Км (7-1)(К-1)В + kB ^ + ka BI 


[G-9(/-2)(k-1)(k-2)(j 1) Gi - 2) 887 ] (6,027) 


54), + 2k, aB U; буд + Ki 95051 + 2k, U;B; (5, + 285, |8, 


B, = Cni G-9(k-), 8 
u—— j+k-3 ]+К-1 


Си = к (К-1)В:7°6,, - Кый, (б, +26,,)B; Bua | В; + (k-1)(k - 2) | 
Dy, = kB)" -kU U8) Bi 8 -(7-117-2)| 

Dp Sk (Î-DBE ; Бд--28410,8 28, 8 4(-10-2)| 
б,-0 pentru iz j si б, =1 (simbolul lui Kronecker) 


Solutia nebanala a sistemului (8) se obtine din conditia ca determinantul 
characteristic al coeficientilor lui a, sa fie nul. 


3.Rezultate numerice 


In primul exemplu, consideram p=2 deci doua bare curbe circulare. In tabelul 1 am 
obtinut primele doua frecvente-parametru pentru diferite valori ale lui m in relatia (1) si 


pentru ky; =k, = k,, =k . In acest caz U, =1, Cj, = Cj, =1. Rezultatele sunt comparate 


cu solutia exacta obtinuta de Sato . 


Numarul 
modulului m k 
10: 104 10° 10° 00 
2 1.7610 1.7911 1.7934 1.7935 
3 1.6321 1.6512 1.6571 1.6572 
4 1.6202 1.6410 1.6413 1.6413 
1 5 1.6201 1.6413 1.6411 1.6412 1.6380 
2 3.7998 4.0517 4.0569 4.0571 
3 3.7972 4.0516 4.0568 4.0572 
4 3.6413 3.9215 3.9231 3.9233 
2 5 3.6410 3.9214 3.9231 3.9232 3.8150 
Tabelul 1: Parametrul О pentru U, =1, C, = С), =1, ky, = Кк, = К, = К 


In tabelul 2 sunt obtinute primele 2 frecvente-parametru pentru U, = 0,168, 
Ch = С, =1, Ки =k, =k, =k. 


Numarul 
modulului m k 
2 0.6625 0.6636 0.6637 0.6537 
3 0.6211 0.6234 0.6235 0.6235 
4 0.6103 0.6193 0.6193 0.6193 
1 5 0.6103 0.6193 0.6193 0.6193 
2 1.6637 1.7213 1.7229 1.7229 
3 1.2135 1.2211 1.2213 1.2214 
4 1.0991 1.1011 1.1021 1.1022 
2 5 1.0817 1.0917 1.1021 1.1022 


=k, =k, =k 


Tabelul 2: Parametrul Q pentru U, =0,168 , Су, = С, =1, ky, А 


Se observa ca viteza de convergenta creste odata cu cresterea numarului m de termini 
din formula (1), fiind ceva mai lenta pentru U, - 0,168. Solutia este insenzitiva la 


cresterea k pentru 10°. In tabelul 3 sunt date primele doua frecvente-parametru pentru 
diferite valori ale Іші U, si comparate cu cele date de Sato. In acest са? valorile sunt 


=k, =k =10° si m=5. 


calculate folosind k,, = ky 


1 


Митаги! 
modulului U 
1 
0.630 
0.437 
0.366 
0.213 
1 0.168 
1 
0.630 
0.437 
0.366 
0.213 
2 0.168 


2 


Prezenta 
lucrare 
1.6337 

1.317 
1.119 
0.871 
0.7131 
0.6122 
3.814 
2.991 
2.392 
1.695 
1.311 
1.0840 


Tabelul З : Parametrul Q pentru m=5 ; С’, = С), =1; kyi 
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VIBRATIILE UNUI SISTEM DE ВАКЕ BERNOULLI - EULER 
AVAND LA CAPAT O MASA ADITIONALA 


1. INTRODUCERE 


In literatura de specialitate exista un mare numar de lucrari care se ocupa cu studiul 
barelor cu o masa aditionala la unul dintre capete care iau in considerare factori complicati 
cum ar fi: forte axiale, arcuri elastice, sectiuni axiale variabile, inertia de rotatie a masei 


aditionale, etc. Astfel Goel [1] foloseste transformata Laplace pentru determinarea primelor 


trei frecvente proprii ale unei bare in consola si o masa concentrata. Metoda aceasta a fost de 


asemenea folosita de Liu [2] pentru obtinerea primelor cinci frecvente proprii pentru bare 


uniforme si neuniforme articulate la ambele capete si avand arce elastice si doua mase 


concentrate. To [3] a obtinut solutia exacta a vibratiilor unei structuri de tip antena, luand in 


considerare faptul ca centrul de greutate al masei aditionale nu coincide cu punctul de fixare. 
Gürgóze si Batan au introdus arce pentru studiul vibratiilor unei bare cu o masa aditionala. 


Kim si Dickinson [7] au analiyat problema barelor Euler avand termeni complicati ca de 


exemplu: mase concentrate, momente de inertie concentrate, suporti intermediari, incarcare 


axiala si sectiune neuniforma Jang si Bert [8] au investigat frecventele proprii ale barelor 


Bernoulli — Euler pentru diferite conditii la limita si doua sectiuni axiale diferite. Laura si 


Gutierrez [5] au obtinut frecventa fundamentala a vibratiilor unor bare de sectiune neuniforma 


folosind metoda Rayleigh — Ritz combinata cu metoda Schmidt. 


2. SOLUTIA EXACTA, CONDITII LA LIMITA 


Ecuatiile diferentiale ale vibratiilor libere a doua bare Bernoulli — Euler pentru 


amplitudini mici se scriu sub forma (fig. 1): 


азу, а?у, 
Eh VA + hA, S = 0 (1) 
E,I, 335 + p A, = = 0 (2) 

272 a= Po 2 Жа 


in саге E, este modulul de elasticitate al barei i (i=1,2), 1; este momentul de inertie geometric 
al sectiunii, р este densitatea de masa, А, este aria sectiunii, L. este lungimea, x este 
coordonata spatiala de-a lungul axei, t este timpul, iar y, este deplasarea transversala a barei i. 


Pentru prima bara incastrata iar a doua avand capatul liber, conditiile la limita se scriu sub 


forma: 


4 dy, (0,8) 
1(0,4)-0, = 
уа(0,4) ах 
709 8) mg 222060 — 
EL + m? та ae ü (3) 
8? y, (Lt) ay, (1,4) а?у (1,4) 
Ez — 5z T (J + md“ 2) SH ЕТЕРІ. + md —— art - 0 


unde L = L, + L. , d este distanta dintre centrul de greutate al masei aditionale si punctul de 
fixare (fig. 1), J = mR? este momentul de inertie al masei atasate m si R este raza de giratie. 


Am presupus ca presiunea de contact dintre cele doua bare este neglijabila. Conditiile 
de continuitate ale deplasarii, a unghiului de panta, a momentului si a fortei taietoare se scriu 


astfel: 


a a 
zz. = = gis (L, t) 


y (Lat) жы y; (L,.t), Jz 


s. B a? 
ЕЁ ag (Lrt) = Е,1, agp (14,4) (4) 
3 3 


a d 
Eih зіл (L, t) = Е, 232220410) 


Definim variabilele adimensionae іп felul urmator: 


(5) 


fh „_ 4 QR  _ т 
ТЕК “ЗО ТЕСР UIT 
Fie 
y (x, t) = Y, (x) cos wt, va (x,t) = Y; (x) сов wt, (6) 


unde Y, si Ү sunt modurile transversale de vibratie iar ¿o frecventa. 

Folosind variabilele adimensionale (5) si relatiile (6), ecuatiile (1) respectiv (2) se 
reduc la relatiile: 
dir 


жад. bir) = 0 9 


4 


ai Y, (x) 


йх^* 


— bY, (x) = 0, 
Conditiile la limita (3) si conditiile de continuitate (4) se scriu astfel: 
Y -аү200, H G) = а"), К) KiYQQ)- К,КЬЗҮД(1)=0 (8) 


Үг (1) + K, (KZ + KZ)biY; (1) - K,K,b3V, (1) = 0 


Solutiile generale ale ecuatiei (7) sunt : 


У. (х) = Cy - cosh b,x + C,,- sinh х + Cy cosb;x + C; -sinb,x, i=1,2 (9) 


3. ECUATIA FRECVENTELOR PROPRII, APLICATII NUMERICE 
Inlocuind conditiile la limita si de continuitate (8) in ecuatiile (9), obtinem un sistem 
de ecuatii in 8 necunoscute: C44, C44, Су, Caz Caz» Caz» Caz Si Саг. Solutia nebanala a acestui 
sistem este posibila daca determinantul sau caracteristic este nul: 
det(a;,)-0, ij=18 (10) 
in care componentele nenule sunt: 


44 = 044 = Aaa = а. = 0; ay =coshyb,; a4, —sinhgb,; аҙ = cosmb, 


аза = sinhmb,; аз; = —coshb,; аҙ; = —5іпһ by; а;- = —cosmb; 
üy = —singb,; a, = bysinh yb; a,, = b,cosh nb, 
аз = —b,singb,; аы = b,cosmb,; а, = —b;sinh пр, 


Пы = —h.coshgb, ay, = —h.singb,.; ag = —b,cos 7b, 


as, = bicoshmb,; a; = bisinhgb,:; as; = —bicosnb, 


az, = —bisingb, а: --аБісовһ ba; az, = —abšsinh nb, 
as; = abjcosyb,; ass = abjsingb,; az, = Бісіпһ gb, 
Gg; = bicoshmb,; аз = bisinmb;; as, = —bicosmb, 
tgs —abisinh mb; а, = —ablcoshgb, a, = —abisin mb 
ашы = abšcos nb;; аз; = (1+K,K,b;)sinh b, — Kb, cosh b; 
ass = (1 + K,K,b;)cosh b, — К.В. sinh by; a;,-(1— K,K; b;)sin b, — K; b; cosh Б, 
ta = (—1 + K,K,bi)cos b. — Kb, sin b, 
ass= (1 — K, K;bz)cosh b, + K (KI + КР) іп b, 
as Í 7 (1 — K,K,bz)sinh b, + K (KĪ + K2)bicosh b, 
а55=—(1+ К,К,Ь2)соѕ b, — K,(K7 + Kz)b3sin b, 
Gag 7 — (1 + K,K;bš)sin b; + K (Kf + K2)b3cos b, 
(11) 

In cele ce urmeaza, vom compara coeficientul b, al pulsatiei fundamentale determinat 
din ecuatia (10) cu cel determinat cu ajutorul metodei aproximative Rayleigh.Schmidt [5]. In 
tabelele 1+4 se face comparatie intre rezultatele obtinute in [5] si rezultatele prezentei lucrari 
pentru diferite valori ale parametrilor Кү, К», Ку si 7). In tabelele 1-3 se considera K, = 0 iar 


in tabelul 4, K, = 0. 


K; 
— —— PVA — — ss -—-—- > s ——  -— > 
0 0.3 0,6 0,9 
Rezultate proprii 1.473 1.414 1.276 1.126 
Lucrarea [5] 2.17 2.00 1.63 1.27 


Tabelul 1: Valorile lui b, (= 9, in [3]) pentru K, = 0. К, = 27/40. n, = 1/2 


K; 
O 03 06 oo 
Rezultate proprii 1.320 1.273 1.136 1.00 
Lucrarea [5] 1.74 1.62 1.29 1.00 


Tabelul 2: Valorile ші b, (= 0, in [5] } pentru K, = 0, K, = 9/8. п. = 1/2 


K; 
O o 06 08 
Rezultate proprii 1.446 1.390 1.254 1.114 
Lucrarea [5] 2.09 1.93 1.57 1.24 


Tabelul 3: Valorile lui b, (= 0, in [3]) pentru K, = 0, K; = 51/80. n, = 1/4 


К, Rezultate proprii Lucrarea [5] 

K. 0.6 0.8 1 0.6 0.8 1 
0.4 0.974 0.927 0.877 0.95 0.86 0.77 
0.6 0.905 0.871 0.836 - 0.76 0.70 
0.8 0.873 0.835 0.793 - - 0.63 


Tabelul 4. Valorile lui b, (= 2, in [5]) pentru K; = 9/8, m, = 1/2 


4. CONCLUZII 


Lucrarea analizeaza solutia exacta a problemei valorilor proprii ale vibratiilor a doua 
bare Bernoulli-Euler cu o masa aditionala. Este data solutia exacta pentru о bara in consola iar 
a doua cu capatul liber, sectiunile celor doua bare fiind diferite. Solutia data este adecvata 


scopului propus. Aproximatiile din bibliografie confirma rezultatele metodei excate. 
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INFLUENTA SARCINII AXIALE DE COMPRESIE ASUPRA FRECVENTELOR 
NATURALE ALE GRINZILOR DE TIP TIMOSHENKO 


In cele ce urmeaza se abordeaza ecuatiile de frecventa si modurile normale de vibratii libere 
de incovoiere alre grinzilor finite, inclusiv efectul de deformare prin forfecare, inertia prin 
rotatie si forta de compresie axiala pentru grinda de tip liber. 


Introducere 
Ecuatiile de miscare pot fi derivate cu ajutorul principiului Hamilton 


BEIGE (1) 


t 


Unde t0 si t1 sunt doi timpi fixate, Ec este energia cinetica a grinzii Timoshenko iar U 
reprezinta energia potentiala. 


Energia cinetica este obtinuta prin formula: 


L L 
2 2 
кү (ране (ata © 
° 2) dt) 2 | dt) 
0 0 
Dar energia potentiala este obtinuta: 
L L L 
2 2 2 
MEME TR EC ms у! Medel) асан (3) 
342 dx J 2 dx J 2 dx ) 
0 0 0 


Unde y(x,t) reprezinta deformarea transversala, уу (x,t) раша de indoire, A zona sectiunii 
transversale, p masa per unitate de lungime, I zona momentului de inertie a sectiunii 
transversale, E modulul de elasticitate, G modulul de rigiditate, k factorul de forma pentru 
sectiunea transversala, N sarcina de compresiune axiala constanta, L lungimea grinzii. 

Cu aceste expressi se obtin, folosind principiul Hamolton (1), urmatoarele ecuatii 
diferentiale partiale: 


2 


2 
eS eaa 9 у) pni ee ш (4) 
dx? dx à? dx 
2 2 
wo [Se Е. (5) 
ах dx J dt 2 


N 
B (6) 


sunt folosite, ecuatiile (4) 51 (5) sunt reduse, respective 


2 
2b пу. АВТ, Di 2:2. 
UST Drew sss а суш) (7) 
r 


" ' 2 
у" —- ү'- ау = 0 


(8) 


Eliminand уу sau y din ecuatiile (8) si (9), urmatoarele doua ecuatii diferentiale complete іп у 
51 yy sunt obtinute 


4 4 25252 2 4 
а а b d d d 
2 (a? H 52) 2-2 гу y оо ЫГ b Уг 0 (9) 
d dx? dt a г dt 2 dx? dt Š 
4 27-12 4 
dw (2,424 v b dw 24 y а2ъ24 \ Lg (10) 
ах ax? dt i. dt dx dt 5 


Pentru a gasi modurile de oscilatie in cazul grinzilor libere, conditii de limita trebuie 
impuse: 


w(0,t)=dy(0,t)/dx , dw(o,t)/dx=0 , w(L,t)2dy(L,t).dx-0 , dy(L,t)/dx=0 (11) 


Solutii 


Procesul de separare a variabilelor poate fi aplicat in (9) si (10). Astfel,vom separa 
variabilele in ipoteza ca: 


yGut)-Xi(x)T(t) , wQGut)-Xo(x)T(t) (12) 
Solutiile pentru (9) si (10) pot fi gasite ca 


X1(x)=Cicos qx + Casin qx +Cacosh sx + Casinh sx (13) 


2-2 2 2.2 2 
XQ) 29-Р - Ч (с, cos(qx) — Cysin(qx) t — С cos(sx) + C, sinh(sx)) (14) 


unde 


2 
iu [22 ñ 82) 7 20 due de Dd ñ pha? " 52) EM (15) 


2 
2472 НЕР ь?) + Pee Е кө? Бе р (16) 


T(t)*p^2*T(t)-0 (17) 
si conditia 


Б>агр (18) 
Ecuatie de frecventa si modurile normale 


Conditiile la limita (11) in ХІ si X2 si ecuatie de frecventa astfel obtinute sunt, respectiv, 
Х2(0)=Х`1 (0), Х`2(0)=0 , Xo(L)-2X" ((L) , X`2(L)=0 (19) 


Е - 2-cos (qL)-cos (sL) + | -sin(qL)-sinh(sL) -0 (20) 


Radacinile еспайе de frecventa Pi і-1,2,3... da valorile proprii ale problemei. Functiile 
proprii corespunzatoare modurilor normale X li si X2i, pot fi obtinute in felul urmator: 


g^ (> - (ча. 2 5 


Xi) - K, cos (83) + Ө "зіп (ах) - : 7 cos (s;x) + талын (six) (21) 


4; 


2 os : 2 s. а? : 5 s. 
Хх, = шыш н 2 _ veu s (5) + se (5х) (22) 


ШЕТТЕП 


| „а oy + ыг Q3 


L 100 200 300 
г^2 
125/4 | 1652.478 | 478.684 | 217.928 
125/6 | 1434.980 | 399.392 | 179.825 
125/8 | 1286.600 | 349.812 | 156.547 

L 100 200 300 


r^2 
125/4 | 3330.838 | 1163.33 | 567.843 

125/6 | 3078.249 | 1005.87 | 478.187 

125/8 | 2873.381 | 899.074 | 420.927 
L 100 200 300 


r^2 
125/4 | 5052.785 | 1972.37 | 1025.31 
125/6 | 4803.226 | 1761.00 | 883.702 
125/8 | 4587.763 | 1606.52 | 788.381 
L 100 200 300 


r^2 
125/4 | 6585.670 | 2813.63 | 1541.02 
125/6 | 6445.301 | 2578.21 | 1358.51 
125/8 | 6264.249 | 2394.56 | 1229.04 


Tabelul 1: Frecventele р1,р2,р3 respectiv p4 pentru grinda de tip Timoshenko si с-/0 


Cu Xli si X21 cunoscute, deformarea totala si panta de indoire pot fi scrise dupa cum 
urmeaza: 


1р; 1+0; 
v6.9 = Y (axe | ) 


t (24) 


w(x,t) = Ys d Pi 21 


t 


Constantele aj, bj, aj, Bj trebuie sa fie evaluate plecand de la conditiile initiale date. 
Exemple numerice 


In multe lucrari se afirma ca prin neglijarea efectului de deformare prin forfecare si a 
inertiei prin rotatie, frecventele naturale astfel obtinute nu difera foarte mult de la frecventele 
date de grinzile Bernoulli-Euler.Dar exista grinzi (structura "D", “Т”, “U”, ...) unde 
frecventele sunt considerabil mai diferite.Pentru a justifica afirmatia de mai sus luam in 
considerare o grinda de otel cu structura “Т”. Presupunem: Е-2,1 [ECUATII SI PUTERI]. 

Primele patru frecvente pentru diferitele valori L si r2 pentru grinzile de tip Timoshenko 
sunt date in tabelul 1 (с-0) si 2 (с-0) si pentru grinzile de tip Bernoulli-Euler in tabelul 
3(c=/0) si 3(с-0) 


L 100 200 300 
г^2 
125/4 | 1656.648 | 483.480 | 222.919 

125/6 | 1438.563 | 403.280 | 183.922 
125/8 | 1289.797 | 353.304 | 160.132 
L 100 200 300 


r^2 
125/4 | 3335.258 | 1169.143 | 574.190 
125/6 | 3082.233 | 1010.857 | 483.514 
125/8 | 2877.044 | 903.505 | 425.608 
L 100 200 300 


r^2 

125/4 | 5056.410 | 1977.884 | 1031.684 

125/6 | 4806.657 | 1765.879 | 889.170 

125/8 | 4590.020 | 1610.945 | 793.248 
L 100 200 300 


r^2 
125/4 | 6588.657 | 2518.688 | 1547.132 
125/6 | 6448.271 | 2582.784 | 1363.857 
125/8 | 6267.117 | 2398.786 | 1233.859 


Tabelul 2: Frecventele p1,p2,p3, si p4 pentru grinzi de tip Timoshenko si с-0 


Concluzii 

Importanta tratamentului prezentat aici rezulta din faptul ca ea se bazeaza direct pe ecuatiile 
care rezulta din cele de mai sus atunci cand yy si y sunt eliminate.Consecinta principala a 
abordarii in aceasta lucrare este ca efectul de deformare prin forfecare si de inertie prin rotatie 
nu pot fi neglijate.Astfel, este in contrast cu mai multe tratamente obisnuite, care se bazeaza 
pe ecuatia de deplasare in cauza obtinuta prin eliminarea tuturor variabilelor cu exceptia y(x,t) 
din setul de mai sus.In cazurile grinzilor de tip Timoshenko si Bernoulli-Euler, efectul fortei 
de compresie poate fi neglijat pentru toate frecventele. 


FRECVENTELE UNEI BARE TIMOSENKO CONTININD UN NUMAR 
ARBITRAR DE MASE CONCENTRATE 


In cele ce urmeaza, ecuatiile de miscare sunt obtinute pentru vibratii libere ale 
unei bare uniforme avand un numar arbitrar de mase concentrate. In vibratiile de 
incovoiere a barei sunt luate in considerare efectele purei deformari si inertia de rotatie. 


1. INTRODUCERE 


Este cunoscut ca in teoria clasica a barelor Bernoulli-Euler nu se tine seama de 
efectele flexibilitatii, a forfecarii transversale si a inertiei de rotatie. Efectul forfecarii si 
a inertiei de rotatie asupra frecventelor proprii nu se poate neglija intotdeauna. 

In aceasta lucrare ecuatiile de miscare sunt scrise cu ajutorul legilor dinamicii, a 
rezitentei materialelor si a elasticitatii. Solutiile ecuatiilor de miscare sunt obtinute 
pentru doua ecuatii diferentiale (in deflectie transversala si unghi de incovoiere). 
Ecuatia frecventelor este transcedentala si dificil de rezolvat. Se de un exemplu numeric 
pentru primele doua polsatii ale unei bare in consola cu o masa concentrate situate la 
diferite distante de capatul incastrat. 


2. Ecuatiile diferentiale ale miscarii si solutii 


Ecuatiile vibratiilor de incovoiere ale unei bare finite, uniforme continand un numar 
de mase concentrate sunt: 


» 2 
шин ми? "| (x- x) 2 =0 (1) 


Ox 2 
i=l ot 
oM 8^vy 2 
Е (2) 
бх st? 
52 
М = -El 
бх (3) 


(4) 


Unde y(x, t) reprezinta deflectia transversala, w(x, t) unghiul de incovoiere, T(x, t) 
forta transversala de forfecare, M(x, t) este momentul de forfecare, E modulul de 
elasticitate, G modulul de rigiditate, I momentul de inertie axial al sectiunii transversale 
de arie A, p masa unitatii de lungime, k un factor numeric al sectiunii transversale, m 
i(i=l,n) o masa actionand іп punctul de abscisa х: e(0, 1), 1 lungimea Багет, б(х) 
distributia lui Dirac. 

Prin eliminarea lui M si T din ecuatiile(1-4), obtinem: 


Aplicam metoda separarii variabilelor pentru ecuatiile (5) si (6), 


presupunand: 
убы) = хде" ylath=x(ye™ i= үлі (7) 
Folosind notatiile: 
2 2 2 
2_ Рр 2 p ‘A р 2 p “р 
a = ; = emm c s = — 8 
KG’ ? EI : po M 


бі а transformarea Laplace ecuatiilor astfel obtinute, deducem: [12]: 
42 


хү(х) = (5? -cosh(ox) + a 2 сов(рх)). x (0) + т. :sinl(ox ) pos SIP) | х(0) + 


e 


(cosh(ax ) — cos(px ))-x,-(0) + = sinh(ax ) — m sin(x) pa (0) + 


b° ЭЭ маа) (2 заво [x — x) - <-ішбі(х- x) Hix- x)- 
T 


Т : 
92 (Мех "E 


г = \ 
- sinha (х- х) - | wa sinp (x- x) j {H(x - xi) |: (9) 


T БЕРЕСІ 2 cosh (ах) + a 2-cos (px ))-x,(0) + 
a + 


2 2 
БОЛН - —— а + (cosh (ох) — cos (px )):x,:(0)+ 


- -sinh(ox ) - — L sin(x) pe (о) + (p2 -ac 2). 2. (Mixx): 


i=1 


(Mi-x,-x;) . 


n 
(=s sinha ( (x = xi) - sins (x- ху) xc) о (м; к xi) i 


-sinha ( (x — xi) _ sine (x - ШЕ - xj) + b*. 


Me 1 


1 
a 


1 1 


cost (x-x) - eom (x- x) у (vex) [sin (xx) sim (i 


(10) 


2 
My ГЭС, po esa же (a2 - 02) + 4b? , 


2 
Do scd e a (42 - 2) + 4b2 (11) 


lar Н(х) este distributia Іш Heaviside. 
Intre coeficientii xk (0), хк`(0), xk”(0), xi (0), К = 1,2 sunt relatiile: 


Р, га b?) x' (0) + ax" (0) 


2 
x0) = ————, » XW*0-a2mx(0-x"(0 s 
b-ac 


x"(0) = a x'(0) +x" (0) x"5:(0) = (v? -a с?) хі-(0)- c^ x", -(0) 
(12) 


3. Conditii la limita pentru bara in consola. Ecuatia in frecvente 


Pentru bara in consola, conditiile la limita sunt: 


хү(0) = xo) х1) = х8) х2(8) = 0 


(13) 
iar relatiile (12) devin: 
4,2 
202) 2 
х", = | : "Х11(0) 5 x" 2 (0) - лет; 
2) А 
х"5(0) = -с х"(0): Х5(0) = x"/(0) (14) 


Solutiile (9) si (10) sunt іп acest caz de forma: 


x(x) = У [U (x) xi (xi) | + У [Usi:CO ха + 


іші 1= 1 


У [Us x (x)] ȘI | Uii-G9 x" (xi) | (15) 


1-1 іші 
n 


х) = У; мо У [Узг@9-х\(җ)] + 


ї=1 ї=1 


У [Уз xo (xi) | + У |Уағ(д-х% (x) | (16) 


1-1 1-1 


Unde funvtiile Ок (х) si Ув(х), (К-1,4,1- 1,n) sunt date. 
Aplicand derivatele de existenta a solutiilor xi(xi), х 1(xi), xz(xi) 81 х”>(хі) 
netriviale, obtinem ecuatia in frecvente: 


» Ux) U»(x)-sg Us(x) | 048) | 
V(x) У) М(х) әу Vail x) 
V"(x) УМ) Мх) Wail x) -ij 


Cu x1(x), x2(x) cunoscute, deflectia totala, unghiul de incovoiere, momentele si 
fortele pot fi determinate. 


4. Exemple numerice 


Consideram o bara de otel de profil “Т”. In acest caz, efectele defOrkatiei de 
forfecare si a inertiei de rotatie asupra frecventelor nu pot fi neglijate. 
Presupunem E= 2,1 x 1076 [N/cm^2], G=3/8 E . p= 7,8x 10^-3 [Kg/cm^3], 
А=30[ст^2], I-625[cm^4], K=1/4,4, 1-200[cm]. Se obtine а^2= 4,2755831561 
x 10^-10 р^2[1/ст^2], b^2 = 1,74838891 x 10^-12 p^2[1/cm^4], с^2= 
3,6424769 x 10^-11 p^2[1/cm2]. In cazul i=n=1 si mi =1, 2, 3[Kg], primele 
doua pulsatii sunt prezentate in tabelele 1 si 2. 


Mi[Kg] 
Х ет] 1 2 3 
50 65,77 65,70 66,00 
100 65,82 65,91 65,94 
150 65,84 65,86 65,88 


Tabelul 1. Pulsatia fundamentala Pi 


Mi[Kg] 
Х ет] 1 2 3 
50 386,27 375,20 336,02 
100 361,35 338,29 321,73 
150 381,82 367,92 356,74 


Tabelul 2. Pulsatia P» 
5. Concluzii 


Pentru alte legaturi ale barei, modul de tratare este analog celui prezentat mai 
sus, conditiile la limita (13) schimbandu-se corespunzator. Din ecuatia (17) se pot 
obtine si celelalte pulsatii . In literature sunt putine cazuri cand se tine seama de 
deformatia de forfecare , inertia de rotatie si mase concentrate, dar studiul facut 
este fie incomplete (ecuatiile lui Love), fie cu metode aproximativa (Rayleigh, 
Galerkin) sau prin metode experimentale. 

Soultiile date in prezenta lucrare sunt obtinute pentru doua ecuatii diferentiale 
complete; se stabilesc relatii intre constantele care apar (12); conditiile la limita 
sunt omogene; se obtin ecuatii in fecvente si modurile normale; se dau exemple 
numerice. 


CALCULUL PULSATIILOR PROPRII CU METODA 
REYLEIGH-RITZ LA O BARA TIMOSENKO IN 
CONSOLA 


1.INTRODUCERE 


In multe lucrari se considera ca fortele taietoare si inertia de rotatie se pot neglija in 
vederea determinarii frecventelor proprii. Dar in unele situatii si in diferite conditii la 
limita aceste doua efecte nu se pot neglija. Problema vibratiilor transversale ale unei bare 
Timosenko se afla in atentia multor cercetatori care considera diferite metode de 
rezolvare. Astfel, Goel foloseste transformata Laplace pentru determinarea primelor trei 
frecvente proprii ale unei bare uniforme in consola legata de un arc elastic. Liu si Huang 
obtin primele cinci frecvente proprii pentru o bara neuniforma continand mase 
concentrate si un arc elastic de rotatie. Kim si Dickson folosesc metoda Rayleigh-Ritz cu 
polinoame ortogonale ca functii admisibile. Hamadan si Jubran folosesc metoda Galerkin 
in studiul frecventelor unei bare Timosenko cu o masa aditionala. Gurgoze si Batan 
considera o bara uniforma in consola prevazuta cu arce elastice de torsiune si translatie cu 
mase aditionale. 

In cele ce urmeaza se foloseste un sistem de patru ecuatii diferentiale in patru 
variabile din care se obtine ecuatia pulsatiilor proprii sub forma unui determinant egal cu 
zero. Deoarece metoda exacta necesita calcule mai complicate se foloseste metoda 
Rayleigh-Ritz pentru determinarea pulsatiilor, ceea ce prezinta avantajul ca ecuatia 
obtinuta este algebrica. 


2.METODA RAYLEIGH-RITZ 


Ecuatiile vibratiilor de incovoiere pentru o bara finita, uniforma, sunt: 


ОТ e? 
go Mt 0) 
e 
M --EL S (3) 
Ox 
2 
Mg V sj (2) 


(4) 


in care y(x,t) reprezinta abaterea transversala, y(x,t) unghiul de incovoiere, T(x,t) forta 
transversala de forfecare, M(x,t) momentul de incovoiere, E modulul de elasticitate, G 
modulul de rigiditate, I momentul de inertie geometric la ariei sectiunii, A aria sectiunii 


transversale, p densitatea de lungime, K un factor numeric de forma al sectiunii 


transversale. 


Pentru a aplica metoda Rayleigh-Ritz, se va calcula energia cinetica si energia 


potentiala: 


И 2) Ox 
E E 212 | хд ali ) x 
S дх 24 Ox 


unde L este lungimea barei. 


in relatiile (5) si (6) presupunem functiile y(x,t) si y(x,t) de forma: 


у(х.) =X, (x)sin pt, v(x.t) -X, (x)sin pt 
Prin schimbarea de variabila si de functie : 


х-Іх, X,=—X, 


ве obtin valori extreme ale relatiilor (5) respective (6): 
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р 
Enc. рахта )ах+?- ЯЕ )d 
1 2 
z (х)ах nee х (х) X, (х) | dx 
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(10) 


Folosind notatiile: 


КСІ2 1 12 
Ег? 


2 
r ? 


L 
A 
si cu ajutorul relatiilor (9) si (10), introducem functionala: 
1 1 | 2 1 1 
F(X,,X,)= КЕ cmax -[x idxenf(x; - X.) ӘЛЕГЕ (12) 
0 0 0 0 
Pentru bara іп consola,tinand seama de conditiile la limita 
Х,(0)-Х,(0)-0 , X,(1)=X,(1) , Х,(1)-0 (13) 


alegem functiile X,(x) si X,(x) ca fiind polinoame de forma: 


ЖЕСЕ (14) 
X, (x) = x(x- PS Вх (15) 


unde А si В, sunt constante. 
Din conditiile (Rayleigh-Ritz) de forma 


F(X 
OF ( 6) 0 | ОЕ (Х.Х) | (16) 


A  — дВ, 


Obtinem un sistem liniar de ordinal 2п+2 in А si B, (i - 0, п): 


M[n(wA, - V;B,)- p'sU;A, |-0 
(17) 
> [w,B, +n(V,A, +U,B,)- p^yU,B, |=0 


j=0 


unde coeficientii W, , О, si У, sunt dati de relatiile: 


1 
№, = xil [(i+3)x шағ 2(i€2)x +(ї+1)х |(j+3)xP!+(j+1)x'|dx 
0 


1 
U. =| (x 22 ЧА —2xi? -х "Уа 
0 


(18) 
V; --[(х” -2х” xx") (јез) х? —2(j+2)x™ +(j+1)x х 
0 
Prin introducerea matricilor: 
[w]- (и) : [v] - (U;) 2 1 1-1 7 
(19) 


(А) (ААА) : (8)-(8.8-8,/ 
sistemul (17) se scrie matricial sub forma: 
W]- p's[U УТ А 
wheel ар, | JE т 
nl] [w]« n[U]- py[U] 8) 
Din conditia ca ecuatia matriciala (20) sa aiba si solutii nenule, obtinem ecuatia 


pulsatiilor proprii: 


det у ud Q1) 


пу] ани" 


Ecuatia (21) este о ecuatie algebrica in p^ de gradul 2n +2 іп care precizia 
pulsatiilor p, , 1=1,2,... este proportionala cu cresterea numarului ne №. 


3.APLICATIE NUMERICA 


Se considera о bara omogena de caracteristi: E = 2,1-10" [N / ст? | ; G=3E/8 ; 
p-7,8.107 [kg / cm? ] ; k=1/4,4 Pentru n=2 pulsatiile p, , p, , p,, p, obtinute din ecuatia 


(21) sunt prezentate in continuare in functie de lungimea L si de raza de giratie r^. 


г? 100 200 300 
125/4 300,273 80,122 36,110 
125/6 251,891 65,972 30,079 
125/8 220,971 57,122 25,991 


Tabel 1:ри]зайа p, 


100 200 300 
125/4 1354,123 446,124 214,183 
125/6 1225,112 379,381 179,225 
125/8 1127,231 336,548 157,149 
Tabel 2: pulsatia p, 
100 200 300 
125/4 3010,184 1154,331 603,415 
125/6 2799,838 1000,381 511,213 
125/8 2634,114 901,478 484,124 
Tabel 3: pulsatia p, 
100 200 300 
125/4 4611,921 1951,811 1103,814 
125/6 4412,871 1748,412 969,412 
125/8 4271,324 1632,389 841,321 


Tabel 4: pulsatia p, 


